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ABSTRACT 


A method of estimating and correcting for the magnetic field of 
a dual spinning spacecraft has been developed by employing an extension 
of the dual magnetometer technique of Ness et al , (I97I). This 
new method is useful for those situations in which a magnetometer 
boom of modest length (7-lCfei) is attached to the spinning part of 
a large spacecraft (80010OO Kg) . The pu2!*pose of using a dual 
spinning spacecraft is to accommodate two types of instruments : 

(1) imaging and similar "pointed" remote sensing systems on the 
stationary platform and (2) fields, jjar tides and other in- situ 
measuring instruments on the spinning portion. Present-day 
imaging systems are well known to exhibit large magnetic moments, 
sometimes displaced from the spacecraft center by a significant 
amount. The new method ass’jnes that the stationary part of the 
spacecraft possesses a magnetic field which is represented by a 
combination of a dipole and a quadrupole field. 



INTRODUCTION 


It is well known that the existence of a spacecraft magnetic 
field can he a predaninant source of magnetaneter measurem^t error 
if not properly taken into account (Nesst 1970 ) » This is especially 
of concern for the class of larger spacecraft constructed with 
mininal magnetic constraints, which invariably produce significant 
fields. In this paper we discuss this serious and continuing problem 
and develop a method of reducing the effects of this contamination 
field. Tills new method is an extension of a sch^ne developed earlier 
which uses simultaneous duel magnetonffiter measuranents (Ness et al . , 
1971 ). It is appHcable to dual spinning spacecraft which possess 
both a "stationary" segment and a spinning segment. 

If the net instantaneous spacecraft field in-flight can be 
approximated near the end of the magnetometer boom as that of a 
spacecraft- centered, tilted dipole, the dual magnetometer method of 
correcting the measured field can be easily employed in the same 
manner as was first successfully used on the Mariner 10 spacecraft 
(described in general by Ness et al . , I 97 I 5 and specifically by 
Lepping et al . , 1975). Briefly, at the end of a 6 meter boom, the 
Mariner 10 spacecraft field varied from 1-47 depending upon the 
orientation of the imaging systan (Ness et al . , 1974). 

The magnetic field of an on-board imaging system is known to 
be appreciable, because it contains invar rods used for lens alignment, 
and invar is a magnetically soft metal capable of being easily 
magnetized by the focusing field or other external fields. The 
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imaging syst«n is often offset from the spacecraft center and this 
can be a significant fraction of the bocm length. Under these 
conditions, it is obvious that the possibility of a relatively 
strong quadrupole moment will exist, as seen it the positions of 
the dual magnetometer sensors. 

Ti:e method proposed here employs a centered dipole- plus- quadrupole 
field model for the stationary part of the spacecraft field and a 
centered rotating dipole model for the spinning portion. The 
principal contribution to the quadrupole moment is expected to come 
from the field of the imaging system, but other much weaker 
contributors also will exist. This complex field of the stationary 
part of the spacecraft is seen by each triaxial magnetometer 
as rotating through 36o^ every spin of the spacecraft, since the 
magnetometer boom is attached to the spinning i»rt. Therefore, 
during one spin period both magnetometers obtain many measurements 
of this stationary field plus the ambient field and the field of 
the spinning spacecraft (which appears stationary at the sensor 
positions). It is possible to accurately estimate the ambient 
field for each measurement if one can subtract from the measured 
field an accurate estimate of the total spacecraft field. 

For fluxgate magnetometers, intrinsic zero offsets can be 
determined in flight continuously by use of mechanical and 
electronic flij^ers. The incremental mgnetic field obtained by 
taking the instantaneous difference between the inner and outer 
nmgnetometer measurements is then a function of only the net 
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spacecraft field for each angle 0 at which a n^asurement is ssade 
(0 being the azimuthal angle measured in the equatorial plsine of 
the sjAcecraft), Below we develop a detailed formulation of an 
extension of the dual magnetometer method, called here for short 
the "dipole- q\Midrupole method", and discuss the results of some 
tests of its performance by applying it to many different spacecraft 
magnetic field configurations. 

FORMULATION OF THE METHOD 

A second order spherical harmonic expansion is assumed to 
approximate well the stationary spacecraft field at the position 
of the dimil magnetometers. [This "fixed" field is seen as sp5.nr.iiig 
by the n^gnetometer sensors.] It is assumed timit the magnetometer 
boom is perpendicular to the spin axis of the spacecraft and that 
the intersection of the boom line with the spin axis is very close 
to the spacecraft’s center of gravity. Therefore, this intersection 
point is a natural choice for the center of our coordinate system 
and likewise, the spin axis is chosen to be the polar axis. Hence, 
measurements are made at 6 = $ 0 ^ at all times, and thus, in terms 
of Schmidt (1934) normalized coefficients, such an expansion will 
simplify to the form: 

« 2(^) [gjcos^ + hj sin^] 

( 1 ) 

+ 3/2(|)^ [-g® + vI(g2Cos2(^ + h28ln24.)] 

" (f) sj vf[g2C08<^ + h2Sin(^J (2) 


k 


and 


- hj co8^} 

+ ^[gjsina^ - h^ coa2(^3, 


( 3 ) 


where Q and 0 are the polar and azimuthal angles, respectively, 
r is the distance from the spacecraft center to the measurement 
position, and a is effectively the radius of the surliest sphere 
containing the spacecraft proper, i»e. excluding appendates. 

The spinning part of the net spacecraft field appears as a 
fixed field in the sensor frame of reference and is assumed to be 
approximated by a centered, tilted dipole as; 

(s) 


B 


B 


r 

(s) 

6 

(s) 


2 (f)' 


^2* 


( 4 ) 

( 5 ) 


(f)^ ij. 


b:'-' - (-)- (6) 

where ^.^(i =1,2,3) are orthogonal components of the dipole moment. 
We now ccmbine the spinning and fixed fields to obtain the total 
spacecraft field 

(7) 


which in cartesian coordinate representation is 

“ 2(f)^ ('4 + 

+ ^ (89 COS20 + hi Sin2<|» - 


( 8 ) 


B^ - (|) (-£3 + gj sln^. - h^ co8(|>) 


+ (^) ^ (§2 sln 24 > - hj cos2<|i). 


(9) 
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and 

3 it 

\ ■*“ (f) ^ (gj + hj sin#), (10) 

. 0 

where ® gj^ following id^itifications were made: 

t; t t t t t 

K = B » K = B„, and s B . 
r x' 6 z 0 y 

Since both the and ^ terms in equation 8 are independent 
of 0 and depend upon r according to different power laws, the 

and ^ coefficients cannot be solved for separately using in-flight 
data. This would jaresait no problem if the spinning dipole moment 
were negligible, because then 0, and gg could be isolated 
and obtained numerically. In general, this will not be the case, 
but fortunately for many cases in practice, the gg term is found 
to be negligible, as shown in the Appendix. Accordingly then, 
the coefficient win be treated as negligible in the following 
formulation. 

Consid^ the fields measured at the itmer (1) and outer (2) 
magnetometer sensors, which are estpressed as: 

^(1) - t^pd) + I'd) + (01) 

and 

t„(2) - f^pd) + 3^(2) + (12) 

-*t 

where is the nmgnetometer ' s intrinsic zero offset, B is the net 
spacecraft field, and B^ is the ambient field. For accurately known 
zero offsets of the sensors, i.e., ~ difference between 

the inner and outer measured fields is then 

- S^d) - 3 j,( 2) - J'd) - t'(2). (13) 
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•Ince the ambient field cancels essactly. Using our "dlpole-quadrupole 
model” field we finally obtain: 

■ 2R2 + gj coB^ -f hj aln^) (14) 

+ (g2 cos2<^ + sin 24 ), 

- RjC-tj + gj siu4 - hj cos#) (15) 

+ R^i^ (g2 sin2# - hj cos2#), 

and 

ABz • R3‘t^ + R^*^ (gj cos# + sin#), (16) 


where ( 
field, and 

and 


i = x,y,z) are the components of the measured difference 


'‘3 ' - (f;)' 


.a V** ,a 


-(f:) 


(17) 


(18) 


are exactly determined by the known positions of the inner and outer 
sensors . 

The nine unknown coefficients in the difference equations (].W, 

15 > 16) can be obtained in several ways. A computationally simple 
method, shown through simulations to be clearly adequate, is to 
evaluate the three difference equations at four symmetrically chosen 
values of 0(O°,9O°,18O°,27O°, for example) which are solved simultaneously 
to produce (3 x 4 =) 12 values of the 9 coefficients. Hence, three 
are determined redun<tently. [One is free to choose other sets 
of four or more angles for determining the coefficients, but this 
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choice provides simpler analytical fonas.] 

Once the coefficients’ values are determine! for a given sjace- 
craft spin period they are used in equations 8,9,10 for the outer 
sensor’s position, r^, to estinate the spacecraft field at that 
position, S^(2). !Aien, equation is used to estimate the 
ambient field, B^, frtmi measur«i«its from the outer aa^etometer, 
^(2), and the accurately known The method, of course, 

depends on the "stationary" spacecraft field rmnaining q\msi- static 
on the time scale of one spin period. Ej^rience ftrom the Mariner 10 
mission shows that this win normally be the case for typical 
spin periods (6-12 sec). 

EXAMPLES OF MET^D; JUPITER- ORBITER-H^BE 198l/l9^ 

To dononstrate the utility of the irethod, various cases of 
realistically simulated JOP 8l/82 spacecraft sm^etic fields have 
been examined. In the examples we compare: 

(1) the true field with 

(2) the field estimated by the "dipole-quadrupole model" as 
outlined above, and with 

(3) the centered dipole field approximation as employed in the 
Mariner 10 mission. 

The latter field (at the outer magnetometer location) is given by 



as described by Ness et al . « 1971; ^ is given by the first part of 
eq. 13 and a = (^^'^^2^^ known spatial "coupling coefficient". 


8 



In lArticular for the outer mgnetc^eter location, r^, the following 
quantities are calculated: 

^l) the magnitude of the true spac^raft field, S^; 

(2) the magnitiule of the difference between the "dipole- 
quadr\^)ole model" estimated field and the true field, and 

(3) the mapiit\;uie of the difference between the dipole estimated 
field and the true field, S^. 

First, we calculate these qaantifcies functions of 0 for two 
different bocaa lengths, 6 and 9 In both cases the true sj»cecraft 
field is simulated by three dipoles whose characteristics and 
locations are given in Table 1 where is a latitude angle, i.e. 
the inclination from the eqmtorial plane. These characteristics 
are realistic, based on: 

a) our general knowledge of possible field strengths of space- 
craft subsysteos in combination, 

b) cur experience with Mariner 10, and 

c) the results of field mppings of the MJB*77 spacecraft which 
have transpired to this date. 

Figure 1 shows the results of these simulation tests. For the 
case of the 6 m boom, the locations of the inner and outer sensors 
from the spinning spacecraft center are r^ = U.7 m and r^ = 7*0 m, 
consistent for a 6 m boom and a spacecraft radius plus strut length 
of approximately 1 m. This case is designated with a "7" in the 
figure. In the case of the 9 m boom, similarly, r^ = 7.0 m aud 
rg = 10.0 mj this is designated with a "10" in the figure. 
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TABLE 1 


Simulated Field Characteristics 


^maglng S^ystem Dipole. (Angle 
estimates based on J.P.L. 
660-22, JOP 81/82 Orbiter 
Description Document.) 


Stationary "axial” dipole 
moment (center of stationary 
moment is distance D above 
spinning center, along spin 
axis) . 


Spinning £entered dipole 
moment. (This point is at 
intersection of spin axis 
and boom line.) 


Displacement vector between 
the spinning center and the 
imaging system moment's center. 
(Estimates based on JOP 
Description Document; see 
above . ) 


- 

lOOC 

3 

gauss cm 

* 

15° 

in 

fixed 

*^IS 

0° 

system 

l^sl - 

3 

150 gauss cm 

% * 

80° 

in fixed 

♦s ■ 

90 

system 

D 

0.5 m 

I3cl ■ 

3 

250 gauss cm 

’fc - 


In spinning 

♦c * 

45° 1 

system 

R » 1. 

6 m 


- 45*^ 



♦ r-0° 
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For the 7*0 m case (thin lines) figure 1 shows an unacceptably 
large (thin dotted) for the entire 360° range of 0 z^aching 
w 1 7 at 0 ° 0^. This must be consid«r^ an error field if a dual 
magnetca^ter method is not employed* Sp (t^n dashed) shows a 
distinct improves^t over for all 0* (thin solid) is < 0.1^ 7 
for all 0 representing an improvement over Sp of better than e factor 
of » 3 or an improvement over of better than 6 at 0 -- 0°, the 
worst case iK)sition. 

Similarly for the 10.0 m case (heavy lines) we see that for 
0=0° tl^ maximum (heavy dotted) is 0,29 r, Sp (heavy, dashed) 
is 0.07 7 and Sq (heavy solid) is 0.^ 7 giving impro'^^aent factors 
of Sj/Sp =4,4 and S^/S^ = l4. So, even in the case of the longer 
boom (standoff 10.0 m) for this spacecraft field, is unacceptable 

(i.e. > 0.025 7) at all 0, and Sp is unacceptable for almost half 
the cycle. Only S^ remains below 0.025 7 for all 0. As can be seen, 
the 0=0° position is where S^ and Sp are maximum and therefore 
is the worst case position; this is true for any simulated field 
since the dlsplac^ient vector of the imging systaa is chv.«en to 
be in the 0=0° meridian plane. Therefore, the reiAtive merits 
of the two methods can be determined satisfactorily by comparing 
their results for tl^ 0=0° position only. 

In Table 2 we show the results of ten tests where the input 
values of and Fg are varied but all the other input parameter 

veG-ues are the same es those used in obtaining Figure 1. Odd and 
even number tests correspond tc short (6 m) and nxjder^te (9 m) length 
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booms 9 resp^tively. first two rofws in tbs tab2«, tests 1 aM 
2 , are identical to those in tl^ figure for all input parameters. 
Notice that the r^ = 10 m cases produce values of < 0.025 7 > 
indicatii^ that a canbination of a nasderate length bocaa and the 
new dipole-qimdrupole method viU alirays yield acceptably snail 
errors. Tests 1 and 3 differ caaly in that is slightly larger 
in test 1 where expectedly is slightly larger, but Sj^ is »iual 
for the two tests and likewise for S^. The same commait holds for 
tests 2 and 4. In fact, ^ and Sq are independent of since 

either method exactly subtracts out the spinning centered dipole. 

In reality, how well this holds depends on how well the centered 
dipole approximates the spinning spacecraft field. 

Notice that Sj^ for test 1 is approximately equal to for 
test 2, and likewise for all even-odd test pairs. This demonstrates 
for the particular simxilated net spacecraft field chosen, that 
applying the centered dii»le field approximation of Ness et al . 

( 1971 ) is equ: valent to lengthening the boom by 3 m and applying 
no correction method. Similarly for test 1 is approximately 
one half for test 2 and likewise for all even-odd pairs. This 
means that applying the new dipole- quadrupole method of approximating 
the field yields a decrease in error by a factor of about 2 over 
lengthffliing the boom by 3 m and applying no correction method. 

The odd numbered tests, where the boom is short, yield modest 
improvement factors S^/S^ of « 6, and the even number tests, where 
the boom is of moderate length, yield large improvement factors 
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Spacecraft Field Error Estimates 
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«15. This shows that the ratio (as well as 

itself increases as the boom is length^ed. Just as the values 
of Sjj and accordingly decrease. This illustrates the value 
of using the longest boom technically and economically possible, 

SUMMARY AHD CONCLUSION 

An extension of the dual nagnetoraeter method of Ness et al« (1971) > 
applicable for data from a dual spinning spacecraft, has been described. 
!Rie schme exploits the facts that: (1) the stationary segment of 

the spacecraft contains a nagnetic field source (imaging system) 
which is strong and significantly offset from the spacecraft center 
and (2) the spinning segment is the platform for the magnetometer 
boom. Hence, dual magnetometer data used in conjunction with a 
dipole- quadrupole net spacecraft field model enables an accurate 
estimate of this field once every spin period; this field is expected 
to be quasi- stationary over a spin period. The estimated spacecraft 
field is then subtracted from the total measured field (spacecraft 
plus ambient) providing accurate estimates of the aiobiait field for 
each detail sample. The method used in association with standard 
instrument zero level correction techniques should reduce the net 
error on the ambient field estimate to » + O.097 per axis when a 
modest length boom (7- 10m) is employed. Likewise the method 
provides the ability to achieve moderately good accuracy even for 
the case of the short boom (6 m) under those conditions when the 
longer boom is technically or economically prohibitive. 
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APPENDIX 


An estimate of the error caused by neglecting the tern containing 
the g^ coefficient In equation 8 Is made here. For completeness, equations 8, 
9, and 10 are examined to ascertain the relative contributions of all 
terms, but Initially only the field due to the imaging system Is con<* 
sldered. 

The spherical harmonic coefficients (g^^, h^^) for an expansion 

made about a center, C, are related to the components (M , M , M ) of 

X y z 

an offset tilted dipole (which represents the imaging syst^) by the 
following transformation (Bartels, 1936) , given here to second order: 

h\-My, 

g ®2 * (2 Az - Ax - My Ay)/a, 

= (M Az + M Ax) /^a, 

^ x z 

g ^2 * ^ - My Ay) /v^/a, 

h ^2 * (My Az + M^ Ax) v5/a, 
and 

h ^2 * (My Ax + Ay) 

where R » ig the displacement vector between the point C 

and the center of the imaging system's moment, and a Is a length scale 
normalization factor. These transformation equations are now used with 
the specific characteristics given in Table 1 for the JOP imaging system 
to estimate the relative order of magnitudes of the terms in equations 8, 

9, and 10. In the following expressions, the brackets represent the 
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relative order-of-magitltudes (i.e. haroM^nlc amplitudes) and the sub- 
scripts denote the dipole terms (D) or quadrupole terms (Q), excluding 
the g® term which la presented separately. First the r 2 ■ 7 m outer 
magnetometer standoff position Is considered. The spacecraft field 
components are then estimated: 
r„ 3 

«*3-[1930]j + [7001 + [lOSlO 

C bM . [966]^ + [«7l 

^ y y 

C(^)^ b']. [259]p + [605 3,^ . 

2 Z 

The g® term Is the smallest and is ^ 1/18 the D term and ^ 1/7 the 

2 X 

0^ term. When r 2 * 10 m we obtain: 

[990 [327] , [424]^ . and [ 76]^0, 

where the dipole terms remain the same as for the r 2 ~ 7 m case. Again 
the term is the smallest and is ^ 1/25 the D term and ^ l/ 7 the 
0^ term. 

In a similar manner, the absolute values of the constant g^ term 
for the two cases can be obtained, and are 0.048Y and O.OllY for r 2 * 7 
and 10 m, respectively. Thus, for the imaging system alone, and for the 
characteristics given in Table 1, the g^ term is clearly negligible. 

For the centered spinning dipole of Table 1, It = 0, which yields a g^ 
term of zero. Also for the stationary "axial" dipole of Table 1 
the resulting g^ term Is considerably smaller than the g^ term of the 
Imaging system, because Mg « M^g, (D») AZg ^Az^g/2 and Axg “ Ayg = 0. 
Hence, neglecting the g® term for the net field of the JOP spacecraft 
Is justified. 
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FIGURE CAPTION 

Figure 1. Simulated spacecraft magnetic field error estimates as a 
function of the (equatorial) azimuthal angle 4>. The imaging 
system is located at (^> » 0°. The designations 7 and 10 refer to 
the net standoff distance of the outer magnetometer for the two 
tests, and S^, S^ and S^ are defined in the text. T i gauss. 
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